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INTRODUCTION 


There is a great deal of current interest in nonlinear 
phenomena. The mathematics involved in describing such 
systems is considerably more complicated than if one is dealing 
with linear systems. This paer is a survey of various methods 
of solving differential equations which arise in connection 
with nonlinear vibrating systems. No exact method of solution 


has been worked out, but many approximate schemes are in 


existence. An ettempt is made to state as clearly as possible. 


the mathematical techniques that are involved, so that the 
reader may readily apply them to other situations, as his 
interests direct. Four of the more common developments are 
considered in this paper. 

Dr. P. O. Pedersen worked out a method whereby a 
Fourier Series is assumed as a solution to the differential 
equation which describes certain nonlinear vibrating systems. 
Since the Fourier series is assumed to be a solution, it 
should give rise to an identity when substituted into the 
differential equation. When one makes use of this condition, 
end assumes that the oscillations are small and the systems 
are not extremely nonlinear, a set of equations are obtained 


which give sufficient data to calculate the coefficients of 
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the first few terms of the series. This gives an approximate 
solution for small oscillations. 

For systems which lead to such complicated equations 
that numerical analysis procedures ere the only practical 
avenue of attack, the method of differences may be used. The 
advantage of this method is that the final result will be a 
function rather than just a number. The dependent variable 
must be calculated for several values of the independent 
variable in the region of the curve which one is working. 
From this data all derivations may be calculated anda 
Taylor's series formed. When this Taylor's series is trans- 
formed into a Fourier series the solution will be in a form 
in which subharmonic resonance may be investigated. 

Another powerful tool in attacking this type of 

problem is the LaPlace transform. This is the first method 
to use initial condition. Here small nonlinearity is assumed 
and if "a" 4s the coefficient of nonlinearity, the solution 
is assumed in series form and developed in powers of "a", 
Each epproximation is complete within itself, and better and 
better results ere attained when higher approximations are 
used. L. A. Pipes uses three approximations for engineering 
work. However as "a" is not very small compared to unity, 
more terms must be used in order to attain the desired 


accuracy. 
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A similar type of solution may be obtained without the 
use of the LePlece transform. The iteration method considers 
@ reasonably generel nonlinear differential equation and 
shows that its solution will contain a subharmonic term. Here 
initial conditions are used. This is not the most general 
type of equation, so that it siahat be assumed that because 
the equation is nonlinear subharmonic resonance will be present, 
However, in all cases that I have considered a subharmonic 
term will be present. 

These methods form a crossection of the more elementary 
developments in this field. Application of these methods 
may be made to many practical vibrating systems, both 
mechanical andelectrical. An example of such a system would 
be found in the design problems involved in the construction 
of loudspeakers. 

A thorough understanding of the ideas and techniques 
involved in this paper will aid the reader in undertaking 


more advanced treatments in this and allied fields. 


Subhermonics in Forced Oscillations in Dissipative Systems 


The following equation affords the solutions to all 
systems, with one degree of freedom. 


rs) 
md % Peek sx = F, sin 2wt 
dt at bE . 


Where we assume 
vad 
= M se 5 ae a a ) 
= 1 (1+ ax + BX" toe. ) 
Beyer Be. (lox $AX-+ 20. ) 
Fos Fo (!)+g)x +EoX© + eel 2 
Assume (1) small emplitudes 
oo 
(2) x = Bo, yan An sin nwt + By cos nwt 
Z 
nel 


((a)) all conditions of convergence are 
met. 


The equation would be, 
Mo 8°x +rp GX 48% = Fo (1+ gx) sin 2wt. 


ate at 
The frequency of the driving force is 2 f = 2w 
f-=zw 
— 
let 
2¥ = ro §o.< 4 “o ho = 81% 
BO Mo ow 
Wo = 80 Po = Xe y= W 


What conditions must be met among the coefficients so 


that 
sin nwt + Bn cos nwt 


9 
= ¢2 Bo oe a An 
& n=l 


is a solution? 


Differentiation may be performed inside the summation. 


Thus Se 
ax = es nw | Bn cos nwt - By, sin nwt | 
at n=l 
2 — oof 
a°’x - -n“w A. sin nwt B cos nwt | 
Pa 
Ne 


Substituting these results into the differential 


equation gives; 
[=e] a 
Bo WG - ya (wg - n°w*) A, - 2¥nw Bal sin nwt 
z Nel 
Ls (we - nw") B, + 2¥nw A, | cos nwt = {Go sin 2 wt 
0 n n 0) 


=) 
Ss Oo ae ee =, sin nwt + B, cos nwt 
UOTE laa gil natlald 


n=l 
It remains to put the right side of the equation into a 
more tractable form. This can be done by using the formulas: 
sinx sinp = 4 | cos (a-f) - cos (“tA )| 


sinx cos@ = 3 [ ean (x-8) + sin («+ B)| 


let 
e=- 2 wt 


eg = nwt 


thus 
A, Sin 2wt sin nwt - An [ coe(2 - n)wt - cos(2 +n)wt| 
rAd 


and 
By sin 2w cos nwt = Bn | sin (2-n) wt +sin (2+n) wt| 
ra 


thus the right side may be written 


eo | 

[*0 ; nema 2wt + Gog, 7. An | cos (2-n) wt - cos(2+n)w 
2 Nel 2 

Bn | sin (2-n) wt + sin(n t+ 2)wt | 

2 


Divide the whole equation through by wh and substitute 
No» E , and Po, for appropriate quantities. 
= Bee 2h 
Bo, 2 [(1 - n*8*)an - 2Po fn Bnpin nwt + (1-n"§ )By + 2nPof An] 
2 


Nel 
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Aa 
cos nwt = E + hobo] sin 2wt + ho ya A, |eos(n-2)wt ~cos(n+2) wt] 


W n=l 


By [ sin (nit 2). = sini (n -2)wt | ; 
Since this is to be an identity demand that the co- 
efficients of like terms be equal. The following system of 


equations results: 
70 - Ho Ay - 0 


(1 - g2) Ay - 2p9 By - ho (By - Bg) 


= 0 
(1 - 8°) By + 2p9 BA, - ho (A, + Ag) = 0 
(1 - 48) Ay - 479 BBy - Ho (2Bp - By) = G0 /uk 
(t S48} Bot 4Dq FE Ay - hy Ay = 0 


=| ae SS 
(l= 98) As.- 6p,€Bs - hy (By - Be). = 0 
(205 oP jong peepgea, tigi (Ay snag) <0 
[1 miGhe)"a, 2 Gopcby = hp (Bo °- Bg) = 0 


t 
> 


(1-168) B, + 8p,§A, + hy (A, - Ag) - 0 
etc. 
Note that in any one equation only even subscripts or 


only odd subscripts appear. 


Thus the above system of equations may be split into 


two sets of equations in the following manner: 


Set I: gibo. 
2Bo - hoAg = 0 
(1 = 48°)AQ - 4rg8Bp - ho (2B - By) = Sosy? 
(1 - 48°)B, + 4pok A, - by Ay a 49 
(1 -168°)A, - apoB, - ho(B, - B,) = 0 


+ 


(1- 16&")3, + Bpo9k Ag ho(Ao - Ag) = 6) 


etc. 

Set II. 
(2 reas eeeect Eh ueetolSu > Ba), |= 
(ee Se )B. tty baad he (Al Ae) | 


oro © oOo 


(1 - #5A, - 6pofB, - hy(B, - B,) = 
(lS ehs. + ep tk + hola Ae) x 
etc. 


One solution which would satisfy the equation would be 
An = 


B, =0 forn=1, 5, 5,... any odd integer. 


O = (,8~ ;@) od - 28908 - pa (ge - 2) 
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Oo (ak fA) ofl + eA sors - ae (“te “ ft) 
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Such a solution however, would not show the subharmonics 


A, and By, which we ere seeking. The problem now has become-- 


Find a solution to set II other than 


An = Bn = O for odd n. 


Ze | Ho | Zq “ad 


5. Py << 1 


Assume dee W = Wo . e E paond BE 


Set I becomes: 
(a) 3Bo = hoAo = O 
(b) -3Ap - 4p9Bp - ho (2Bg - By) 


CO/wé 
(c) -5Bo + 4D940 - hoA, = O 


(d) 154A, + 8p9B4 + ho (Bo - Bg) = O 


(e) -15B4 + 8poA4 + ho (Ap - Ag) = 0 
etc. 

Set II becomes: 
(f) 2D9B, + lp (B, - Bz) = 9 
(g) 2ppA, - Ho(A, + Ag) == 
(nh) 8A, + 6p9Bz + ho(B} - Bs) = 0 
(i) -8B, + 6poAz + Hy (A,- Ag) =4 © 


The first few terms of the series should give a 
reasonable approximation. 


bie Pee 8 «el Bo; Al, Bi, Ao, Bo 


Using (a): Ap = By 
2h 


Using (b): Ay 


2 +[% + 429By + oP] where hoB, = 0 
WO 


Using (ce): Ap = SBo Pere ieee 


4P | 
Put (a) and (ce) into (b) 
3Bp = -4{% 4 4P0Be 4 2 ie (2) 
409 we 4Do 
Solving for Bo: 
[ase + Pee] Tt “SR 


From previous assumptions 


4D0 Fs ons = 0 


o @ Bo tI a 4p 9G 


“one 


But Ap = SBo From equation (c) 
4D, | 
‘aan Ap = hed G, 
Swe 
and By = 2hyAp = - ZhpGo 


S 
Sw“o 
For a solution other than zero to exist, the determinate 


of the coefficients of the equations in Set II must vanish. 


6 Ay + (2p0 ho) By + O Ag - hoBsz = O 
(2p9- ho) A, + OB, -hA, + OB, = O 
OA, + ho By + 8Ag + SpoB3 - O 
ho Ay + 6) By + Sp9 Ag - 8B - O 
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fo 
O 2p9 hg fe) -ho 
Roe 29 - Ho O hy fe) a 
0 No 8 SD5 
Ho 6) Po -3 
Expanding: 
(209 Bo) |2p95- hy -hy 0 229 - hy -hy 0 
) 8 6Po| + Ho fe) 8 Ho 
= Oo 
6p ) e9 
(2p9 hho) (2p9- ho) 6p -g | thol2p9 - ho) hy -8 
J né 2P9 ~ Ro = =H Pg 
5 SPo 


thus 
Oe He 7 = cd rd - 
(296 v6 64 ser) | + (2pgh+ng§ )[ - Soho] 
- ho [ 12p9 - 6hgpg + ho |= 0 
arrange this polynomial in powers of Ao 


4 3 2 2 Brio 
-h) + (6p5 = 6P9)h5 + (64 + S6p6 - 12p6 -12p9) ho 


+[O] hp - 28605 - 144n = 0 


or 


§ 
oO 


fe) 


Ae es ae 2 4 
by (64 + 12p6 )né + 256p% + 1445 


let u - h@ 


i 
to) 


uv - (64 + }2p2) u + 256p2+ 144p4 


ree er oe “kee der 
° aL 6 + om - aia 


Pall ‘oy ane, »e 


u = [64 +12p2 | + /[e4 + lep2|* -4(256pF + 144p¢ ) 
ee 
U = (324 6ps] + [2024 + 12875 - 10806 | 2 


Expand the square root by the binomial theoren. 
let c = 128p5 - 108 pg 


al 
1024 +cl® ~ seoxh( dy )e - 32, 6 
pigs ere 
= 52, 128p6 - 108p9 
64 
= + 2p - 4 
- 382 2P6 &7 Po 
U = S2 + 6p - [82 +20§ - 22 v6 | 
16 
ee 
Us he 
: 2 ey 2 
otis 5 | + 2vp 
6) 0 3 a 
or ho = Po [4 + 27 °6] 
. “16 
Expand by the binomial theorem: 
5 
=Ro = =Po + 2 Po 
or, -by.cs 1 teat pf 
— 128 
2P0 
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Using the assumption that Po RL 


Bo = Bifo = 21 
EPo 2Wolo 
or 
Hog, cs 2Po = 2p08,F 
LWoro 
end 
SiFo = 2wWoro 
Evaluating the coefficients for ho = Do: 
It has been previously shown that 
Bo = - 2hoGo 
3 we 
S 2 
Bo = - 260 2p9g1Gomo |= - 46g Pogi 
Mo 
but r= 2 = 225% 
Mo 
re 
Bo = - =) - By = - 4¢,46 
6wo 2Wo 12w6 
Bo = - 8105 
3wG 
Using equation (g) 
20541 - HAL = 0 assuming Az =O 


this leaves A, arbitrary. 


ator = of al ‘$ 
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Using equation (f) 


£P>B, +h, BL - O B., = 0 
(2p9+h) ) B= 0 
Age By — 8) ‘ 
Bo was previously calculated as : | 
Owe 
since pp = X 
Wo 
Bo = -_4¥G 
owe 
Ag was previously calculated as: 
- G 
fe) 
ees 
OwWG 
The solution may now be written: 
xX = Bo 4 Aysin wt + Ap sin 2wt + Bo cos 2wt 
a. 
X = - g, cn 4, 5in wt - Go Sin 2wt - 4¥ Go cos 2wt, 


Consider the apa 
Ho Kl + fix) a a*x + ro (1+ Bx) ak 4 So(1l+ «x)x = F sin out 
B at P at 


Or, using previous notation: 


(1+ px) deg +) 28 (it px)ax + we (1+ «x)x = Gosin 2wt 
at 
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If you plece a further restriction upon y, f ana OC 


such that ia ci then a transformation may be made to an 


equation whose solution is known. This restriction physically | 


says that the coefficient of the linear variation of mass, 
resistance, and stiffness is the same. It may not be a natural 
restriction but it looks as though it could be realized in an 
experimental model. 


(1+ px) d&x y 2¥ Gx + vex] - Gpsin 2wt 
2 ate at | © ° 


d&x + oy Ox + W5X = Go sin ewt 
ate dt Lit Mx 


buteit powers of ve higher than unity are neglected: 


Lon. St] 2 ex 
L+ ux p 


let ya = 8) 


then a“x 2y.0x 4. wex = Go (1+ ¢)x)sin 2wt 
ate at 


This equation has previously been solved in detail. 


Consider the equation: 
My (1+ pxgy + Toll + px)dx + 8o(l+ax)x = Fo(l+g)x)sin 2w 


divide through by mj and make substitutions for Yo. 80 and Fy 
Bo’ By By 


res 2 = 
(1+ pag 2¥(1+ ex) ax + wWo(l+ «X)X = Fo(1 + g)x)sin 2wt 


t 


2) Bria A A, soem nottetteeg4 conde & ees aout 

ns ot sham ed ese NOL sinuO? areas a nerd = bey teat dsua 
vilsoiaydd achiotut ses, atid?  cowodat ats aotiglor’ enon ‘noltteupe 
20am To aottetisv waeall edt Yo das toLtteos ont fins eYyse 


fatstsa s ed ton Yeu ol, ..6mee off. ef seenttite fie (SonsteLeo% 


ng ot bestises ed binds) ot dguods ae adool ai Tud noitolitaes 
| -Lebor fetnemtteqxe 
_twS atag® = bes > ee os eB) (A + f) 
twS fla 08 = xew 4 aie - ze 
xed TS 
:betoslpen ona. bp neat ‘tenyte 43 to lane. tt dud 
xK~ ft = 


a 
two ate(xpa +t) gd S xew y xbye ri x* nontd 
eT Mogi, Sh 


tested at bevios nosd yLavatvery esa noltaupe elt 


wv nealx ps an s siamo gine Epewen ey 
= | e* 7 so tot snotsus tteda oxen an vd conse abtvte 


c oy i 
H a 7 ‘ - i 


tw ve. aa = xh + “2 + ‘ 


SK se | ee ——— i Se ae 


nee le el eee 


Ik, 


a | ED 2 


To a first approximation this may be written as 


e™ a’x , a¥ebX ax 4 wee%x =» Go(l+g)x)sin 2wt | 


ate at 
) 

If we assume = A = 

eh*| atx + 20x we ad = Go (1 + gx) sin 2wt 
ate at 

or 
a + 2¥ ax + Wex - Ge P* (1 + gx)sin 2 wt 
ate dt 0 - 


ox 


let 
2 2 
d 2N a Wo| X 
5 +a + 
E dt ; 


Expand e~P* and neglect terms of 2nd order 


Ox 


Gp (1 - BX) (1+ g)x) sin 2wt. 


- g,[1+ (1 -6 )x | sin 2wt. 
let gi “Be a 
(Jx = a, [1+ ax] Sin 2wt 


This problem has been solved previously. 


The generalized equation is: 
mo (1 +yx) de X 4 Foll+ / px)ax + 8o(1+ <x) = Fo thseyeieieed 
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n=1 k=1 


= oO 
+W6 ie ll Bo + D_ Ansin nwt + Bycos nwt | 
P34 2 Nel 


(~"e) 3d 

2 

‘aca z Asin nwt + B cos sai 0+ a A,sin kwt + B, cos wt] 
N= Ks 
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cO 90 

ps Lure 2 Slue 

ap = n°w°s A sin nwt sin kwt +(n2 + k2)w AB, 
kel Nel 


sin nwt cos kwt + n@w°B B, cos nwt cos xwt | | 


This is the complete expression for the product of the 
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SPECIAL CASES 
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) A= xX = g, = O 
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The solution for these special cases may be obtained by 


using a procedure similar to the problem previously considered. 
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This set of three equations may be solved, when f(x, x) 


is given, for 60) (0), and xl), 
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Since x and x are periodic functions it is possible 
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x = Zi w<(1-n“) * = 


n=2 
90 
+ bncos nwt 
mie w(i-n7) 


Thus it has been shown that subharmonic resonance will 


be present in a large class of vibrating systems. 
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Solution by the Method of Differences: 
2 ra 
° 


e 


fae een ae OK, xt) 


ate 

By the use of numerical methods values of xX may be 
found for given values of t. Using difference methods 
derivatives of x may be calculated to any order. 

Using the above information e Taylor Series may be 
written for f(x). 

f(x) = fla) + (x-a) £ (a) + (ea) 2" (a)+ ... 


let h - x-a 
re) 
(x=) ate Ga) 
may 


n=0 

The calculated values should be taken as close to the 
portion of the curve being considered as possible. 

Many times error is introduced when calculating 
derivatives by the method of differences so that care should 
be teken to check on this error. 

Since a periodic solution is desired we may write: 


a0 
f(x) = Bo, >. Ansin nwt + B cos nwt 
2 Nel 


The coefficients of such a series would be given by the 


ordinary methods. 


Te f(x) «= o n” an) (a) 
n! 


n=O 


then 
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sanlel sin nwt dt = A, f san? nwt at 
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Me . 


= 
ne (m) (2) cos nwt dt = B cos<nwt dt 
mse 


n@p(m) (2) 
m! 


dt = Bo 
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EXPERIMENT I 


ifaf | 


A - Oscilloscope 
B - Neon Bulb 


This experiment is described in Nature Magazine of 


September 10, 1927. It was claimed that subharmonics up to 


1/40 of the driving frequency were observed. 


set up it yielded no subhermonics of any sort. 


When this was 
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EXPERIMENT II 


To Demonstrate the Actual 
Existence of Subharmonic Resonance 


' OO 00 


High frequency generator 

Iron cored inductance 

Air cored transformer 

Variable parallel resonant circuit 


Amplifier and Oscilloscope 


This experiment is designed to generate and detect 
subharmonic frequencies. 

The varying permeability of the iron cored inductance 
produces non-linearity in the primary circuit. This non- 
linear condition implies that subharmonic frequencies of 
oscillation mey be present. 

The tank circuit may be tuned to select any subharmonic 
frequency. This frequency is amplified, put on the screen 
of the oscilloscope, and observed in the form of visible 


oscillations. 


J 


Mathematical Analysis of the Experiments 


Reference Figure Experiment I: 
Assume that the voltage drop across the neon bulb is a 
function of q. 
Kirchhoff's law gives the following: 
Esinnwt - f(q)-g = 0 
g<=+ -flqa)+E aa 
) 


let F(q) = - f@ 


differentiating this equation gives 
l dq = Ma) dq 4nw E cos nwt 


c dt dt 


or (7) 
2 _F (a) | aq - nw E cos nwt 


Qs 
ct 


since dq zi 


at 
i <= nwE cos nwt 1 | 
¢C l=-c FI\q) 
WW) 2) 
4 = nwE cos nwt [2 + eF (a) + of Gaeta | 
Cc 


assume c << 1 
Then the series may be shown to converge and is a 
function of q. Designate this function as Q(q). It can be 


expanded into a Fourier series. 


co 
Q(q) = Bo+ 2 Ap,sin nwt + Bycos nwt 
N=l 


at 


This series is multiplied by cos nwt. 


It 48 known that 


|e 


cos nwt -cos(n-1)wt - 4 cos wt + 4 cos(2n-1)wt 


This shows that subharmonic terms will be present. 
Variation of Experiment I: 
Place an inductance in series with the condenser. 
The equation describing the circuit becomes: 
Esinnwt - f(q,4)-9q -Ldi = 0 
fs) dt 


Here the function describing the a drop across the 
neon bulb may be made more general: ff (q,4 


or Ldi,q = -f(q,q) + E sin nwt 
dt c 
let F(q,q) - -f(q,4) 
Ldgi,qd = F(q,4)+ E sin nwt 
dt .¢ 
but di -< aq 
at ate 
a ae a*a,, 2d sek F(q,q), E sin nwt 
at le L L 


This is exactly the same form that was considered under 


the iterative procedure. Here again theory predicts that 


subharmonic resonance will be observed. 


— = 


Reference Figure Experiment II: 

Assume that the effect of the secondary circuit on the 
> primary circuit is negligible. | 
In general for an iron core inductance i = f(@) 


In this case assume 


f=¢ 4 1¢° | 
L b 


N1+No is a constant N 


Esinnwt =N ad , 441k | 
at Cc 
nwE cos nwt = N ag 1 dq ,R ai 
datz ‘cat. at 
let nwEe = A ; i jeans ey; B =e D 
N Ne N 


A cos nwt - a°g_+Bi +D di 
dt’ at 


Substituting for i gives: 
d&g + D/L + SO ag 4, B/G, BG] = Acos nwt 
at L b dt L ob 


This equation may be arranged in the following manner: 


ap, Eg = - DI 1,30%}48 _ BOO 4 Acos nwt 
dt L Ep: at bp 


This is the same form as was considered in the iterative 


method. This method showed that subharmonic resonance would 


q be present. 
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A Short Introduction to the LaPlace Transform and its 
Properties. 


The LaPlace transform has found much favor in recent 
years as a powerful tool in the solution of differential 
equations. The uniqueness of this tool displays itself in a 
twofold manner. In the first place, only relatively simple 
integration is required to solve a differential equation by 
this method. This is often a great saving of time and effort. 
The second advantage is seen in the fact that the initial con- | 
ditions are used directly in the solution of the differential 


equation rather than considered seperately as is the case in 


conventional methods of solution. 


DEFINITION: LaPlace transform of f(t) is ferPte(siat. 

Thus the transform is a function of p over the range 
of p for which the transform exists. 

There are many conditions to be placed on the function 
f(t) to ensure the existence of the integral,but all of these 
conditions are satisfied by the physical phenomena treated in 
this paper. Thus in all the cases considered here the integral. 
will be defined over the positive range of p. 


Properties that will be used in this paper: 


Let L{t} = LaPlace transform of f(t). 

1. L f{at(t)+bg(t)} = aL {e(t)}ed Lég(t) 

2. L {ah cto). ph L{t(t)} - p®-2 £(0) - p22 a £(0) 
at 


OOO = n-1 
= nit? : 


Examples of LaPlace transforms: 


L {1} = fa e~Ptat - -e7P*| - 1 pro 
° Pp > ~«~@?~P 
~ oo 
L{ eat} . feat et at = fer(r-alat 22 pra 
° ° p-a 


Example of a solution to a simple differential equation with 


given initial conditions: 


a2 = a& y where "a" 4s a constant 


az y - ak y = O 
at 
Initial Conditions: at ts 0 


Jy 
dy = 0 
at 


ot? a pt ae 2 ae 2x0 
e~ - av y = 
[seo Py] 


2) 


oo oo 
e~pt ab gat - of [oP y at = 0 
if dt o 


Integrating the first integral by parts twice 


enPt a2 y at = e-Pt ay | yp enPty| 4 pe fort y at, 
[or ' 


re) 
notice that: 


or 


1. e-Pt ay la - -dy(0O) =0 due to the initial 
at '. at conditions. 
2. pe~Pt y| = -py(0) = =p 
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The equation may now be written in a factored form: 
ad 
(pe - a&) [ev ydt-p20O0 
6 
ers 
ao 
fevt y dt s 
° pe - ak 


This equation says that the LaPlace transform of y - 
p 


Thus by consulting a table of LaPlace transforms you may 
readily determine what function has a transform js 

Ce ae 
This function is the solution of the given differential 
equation: 


y =» cosh at 


Check: 


ad ( cosh at) = ‘a sinh at 
at 


_a@. ( a sinh at) = a” cosh at 
at 


Substituting this into the differentiel equation it is ob- 
served that the equation is identically satisfied. 


Are the initial conditions satisfied? 


1. cosh (0) = 1 sg y(0) 


2. asinh (0) =» © = ay(d) 
at 


Thus this solution satisfies both the differential 


equation and the initial conditions. 
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Application of the LaPlace Transform (5) 
m aex F(x) = @) 
7) As 


zw Kx + bx® 


thus d&x 4 kx, b x= 0 
Hd m m 


at 
where ke w andb= x 
m m 
the equation would be: 
akx 4 wex 4 ~ x9 - O. 
ate 
Initial Conditions: 
Buppose, at t=0 X= 8 
ax = 0 
at 


Teking the Transform of the given equation 


t | 33 + W°x +x | dt =O 
ate 


Using betes eae of the transform as previously 
stated, 
L | aa! = p*L{x} - p flo) - 4 £(0) 
ate at 
x(©) .= .£(0) thus 
p f(0) = ap 
af(o)= 0 
at 
tes E = L{x} 
L { wx} = wx 
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The expression of the differential equation may now be 
put in the form: 


p°R ~pa¢wk+ aL {xl = 0 | 


or 
(p© w-) X # pa- axe 
Let aa Xq + &X] + av Xo + LXy +... | 


Ze 3 


= 
al 


© = Xp + AF, +4 Xp t oC Fy Ts se || 


Here all powers of x higher than the Srd shall be omitted. 
See (De 4 We + %C,+ a0, +405) ( Xo + «X) + x X, + < X) = 
pa - XL { (xp 4% xX, + Xo ae axa)? } 
Performing the indicated operations on the above equation 
and grouping like powers of x gives: 
Kop? + wk + x [eX + pe + weX, + L {xs} | 
aL af pz, + Wok, + CoXp +C)X, + 3 L {x, x6 }] 


+o] px, ¢ we Ext CX, + Cok) + Cok + L{3x)x% + 3x,x8}| — 


If this is to be an identity the coefficients of like 


powers of x may be equated on each side of the equation. 


The first equation would then be 
Kop’ + wok = pa 


There are two things that should be noted about this 


equation: 
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1. It completely determines Xo which is the first 


approximation. 
2. Further refinement of the solution will not change 


Xo- Thus with this method the complete solution is not re- 


quired in order to determine whether sub harmonics are present. 


ice 


From a table of transforms, 


e e 


Xo = @ COS Wot 


This result carries the method asfar as would be 


necessary to determine whether or not a sub harmonic was 
present. 


Application of LaPlece Transform to a Forced System 


2 2 5 
QAX +wx+xx ws A sin 2wot 
dt© ° 


Here the notation is consistent with the last section. 


[= 2) co 
[er Ee + w°X + «x | at. =) on [ePtetn owt dat 


at A 


The left side of the equation will be the same as it was 


in the previous section. 


ao 


[ert a@x 4 w°x + xx] dt = (p® + w°) X - pa +aL{x®} 
i ate 
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L { A sin 2ut | = 2A wo 
Bee ane 

If we expand x, w, and x in series exactly as was done 
in the last section nothing will change except L A sin 2wt 
will be added to the right side of the equation. 

Since only a first approximation is desired to determine 
whether or not a sub harmonic is present only the first term 
of the left side will be reproduced here. 


ou Vomad 
De Xn. Woe. es pa + A Wo 
0 00 | p*+ 4w \ 
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The w will contain w,terms and as such will be necessary 
to obtain higher approximations but at least we will be 


assured of a term whose LaPlace transform is 
pa - L { a cos Wot } 
Pp + W6 
Thus it has been shown that a subharmonic will be poresent 


in the solution of such an equation. 
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Application of the LaPlace Transform to the Case of an 
Unsymmetric Restoring Force. 


m a°y , kx + ax© = Asin 2wot 
at 
ne ay ee ae = Asin 2wot 
ate m m m 
let kX s we a = & A = F 
m m m 
then a@x + wex + x” - F sin 2wit 


Initial Conditions: 


at t = O b 
ax 
at 


a 
6) 


[ort [ ay vi vac?) a : L{F sin 2vot | 


u{ are, 
L { w2x | 
L(x?) 


The equation 


p°x + wx 


(p©+w’) X 


“<< 1] 


Assume 


= p°x - pa 
= wex 
= XL { x*} 


may now be written in the form. 
- patexL {x*} - L {F sin 2vot } 


= patlL { F sin 2wot } + x L{x?} 
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Substituting these relationships into the equation: 


ewe x +t Ki tak 4 KX, + 0k 
PY + Wot XX + XX, x5][ 9 +X X,+ x X,y X. 


pa + L{F sin ewot }+ x L{x?} ° 
Since we are not necessarily interested in a complete 
solution of the equation but only in whether or not a sub- 
harmonic term is present in the solution we shall take only 
a first approximation. 


(po + we) ¥ pa+L {F sin 2w,t } 


= pa L{F sin 2wot} 
- Pear. Wo 2 P~ + WO 


Considering only the first term, 


= a. 
terre 
fe XO = 


& COS Wot + 7 L {F sin uot} 
p+ Wo 


This shows that the subharmonic frequency Wo will be 


present in the solution. 
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ABSTRACT OF THESIS | 


What is subharmonic resonance? Where and why should it 
occur? These are some of the questions which are considered 


in this paper. 


Subharmonic resonance is a phenomenonthat occurs in many 


vibrating systems. It has been observed in fields as widely 
separated as acoustics and aerodynamics. F.R.W. Strafford 
observed second subharmonics in a moving coil loudspeaker, 
while others have observed airplane surfaces vibrating at a 


frequency which was a subharmonic of the engine frequency. 


Essentially, subharmonic resonance occurs whenever a system 
vibrates with a frequency Z times the driving frequency 
where n may be any positive integer. 

At present there seems to be no clear cut explanation 


of this phenomenon. It has been suggested that higher 


harmonics occur in such profusion that among these harmonics 


there may be contained a series, which has, as its fundamental,,. 


a frequency lower than the driving frequency. Due to the non- 


linearity of the system this "fundamental" frequency may be 


excited. There is some basis for this since it is well known 
9 that the fundamental may be filtered out of a sound and that 
the ear, due to non-linearity, will fill in the fundamental. 


Subharmonic resonance does not always arise because of 
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non-linearity. Take the case of a string, one end of which is 


fixed, and the other end attached to a tuning fork. The 
tuning fork may then be set into vibration and will drive the 
string. Since the string can undergo tension only, it is 
excited by the tuning fork only half the time. Thus the 
string will vibrate at half the frequency of the tuning fork. 


The experiment which was adopted to establish the ex- 


istence of subharmonic resonance gave absolutely no indication 
of this phenomenon. It has been suggested that one possible 
reason for the experiment's failure was that not enough power 
was put into the system to excite the subharmonic frequency 

to an observable level. 

The second experiment was designed to eliminate this 
condition. The subharmonic frequency is generated by a non- 
linear inductance, is selected, 1f present, by a tank circuit, 
and is amplified until observable on the oscilloscope. Due 
to the lack of time this experiment was not performed. 

The theoretical aspect of this paper deals with four 
different methods of solving the differential equations which 
arise in the consideration of certain non-linear vibrating 
systems. 

The method which is developed to the greatest extent in 
this paper is due to Dr. P. O. Pedersen. He assumed small 


oscillations and small non-linear terms. The method then 


boils down to substituting a Fourier series into the given 
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differential equation and making certain simplifications based 
on the assumptions. The relationships which are then found | 
enable him to arrive at an approximate solution. The most 
general form of the differential equation is considered first, | 
and then it is specialized and simplified as convenience and 
interest dictate. It was observed that two rather difficult 
cases could be transformed into a simpler case, already con- 
sidered, if the constants in the linear variation of mass, 
resistance, and stiffness were considered equal. This would 
not be expected in the most general case, but it is probable 
that such a condition could be produced in the laboratory. 

The methods developed in the treatment of non-linear 
vibrating systems by using the LaPlace Transform are due to 
L.A. Pipes and may be found his book, "Applied Mathematics for 
Engineers and Physicists." This method assumes that only a 
small degree of non-linearity is present in the system. 

If "a" is a measure of the non-linearity, then the solution 
to the equation is developed in powers of "a". This can be 
seen to converge if "a" 4s small compared to unity . To 
obtain higher approximations using this method one has only 
to add more terms rather than making an essentiel change in 
the method. 

The numerical method makes use of the fact that the 


derivative may be calculated when only numerical values of 


a function are given. Once these derivatives are known a 
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Taylor series may be formed. This is a solution to the 


equation, but it is not the desired form. From the Taylor 
series a Fourier series may be computed. This method would 
not be very prectical unless the equation would yield to no 
other treatment. 

The method used by T. von Karman employs en iterative 
procedure and also assumes a small non-linearity in the system. 
The first approximation is found to be periodic and is then 
immediately expanded into a Fourier series. This procedure 
is somewhat more general than others considered since an 
arbitrary function is maintained throughout. 

It is interesting to notice that in every case considered | 


@ subharmonic term was found in the solution. This fact leads 


me to believe that subharmonics are present in all non-linear 


vibrating systems. 
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